Off-diagonal hyperfine interaction between the 6pi/2 and 6^3/2 levels in 133 Cs 
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The off-diagonal hyperfine interaction between the 6pi/2 and 6P3/2 states in 133 Cs is evaluated in 
third-order MBPT giving 37.3 Hz and 48.3 Hz, respectively, for second-order energies of the 6p 3 / 2 
F — 3 and F = 4 levels. This result is a factor of 10 smaller than one obtained from an uncorrected 
first-order Dirac-Hartree-Fock calculation and used in the analysis of a recent high-precision (< 2 
kHz) measurement of the 6P3/2 hyperfine structure [Gerginov et al. Phys. Rev. Lett. 91, 72301 
(2003)]. The factor of 10 difference has negligible effect on the conclusions of the recent experiment 
but will become important for experiments carried out at a precision of better than 1 kHz. 

PACS numbers: 21.10.Ky,27.60.+j,31.15.Md,31.25.Jf,32.10.Fn 

I. INTRODUCTION II. PERTURBATION EXPANSION 



In the recent study of the hyperfine structure of the 
6p 3 / 2 state of 133 Cs by Gerginov et al. 0, intervals be- 
tween hyperfine levels were measured to an accuracy 
of < 2 kHz, which was sufficient to give, for the first 
time, a non-zero value for the c hyperfine constant. The 
value of the nuclear octupole moment of 133 Cs obtained 
from c was il — 0.82(10) & x fi N , which is about a 
factor of 40 larger than nuclear shell-model prediction 
^n.s.m = 0.022 & x /ijy, motivating a re-examination of 
corrections to the hyperfine constants. One such correc- 
tion is the second-order hyperfine interaction between the 
6p 3 /2 and 6p±/2 states. An estimate of this correction, 
based on an independent particle model of the cesium 
atom, was used in 1]. In the present work, we carry 
out a detailed third-order MBPT calculation and obtain 
corrections to the 6p 3 / 2 hyperfine levels that are a factor 
of 10 smaller than the values used in p|. Revised val- 
ues of the hyperfine constants a, b, and c, obtained using 
the present results for the second-order hyperfine ener- 
gies agree with those reported in Q to within the error 
estimates. However, for future experiments, especially 
experiments aimed at a precision of better than 1 kHz, 
it will be important to use the correlated values of the 
corrections presented here, rather than the larger values 
given in 0. 
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We write the hyperfine interaction in the form 
H u = £(-l)V fc A >Mf\ 

fcA 

(fc) 

where T_ A is an irreducible tensor operator acting in the 

ffc) 

electron sector and is an irreducible tensor operator 
acting in the nuclear sector. The first-order hyperfine 
correction to the energy of a state |1) is 

W [ p ] = (l|# h f|l) 
= E(- 1 ) /+J+F {/ J I ){j\\T^\\j){l\\M [k) \\l)- 

k ^ > 

The nuclear matrix elements are given in terms of con- 
ventional nuclear moments through 

(lI\M^\Ij) = /i 

(ii\mP\ii) = \q 

(|//|A/ (3) |^) = -ft- 

Here, /j, is the nuclear magnetic dipole moment, Q is the 
nuclear electric quadrupole moment, and SI is the nuclear 
magnetic octupole moment. With these definitions, we 
introduce the conventional hyperfine constants a, b, and 
c through the relations: 

a =j 7j {jJ\t£ ) \Jj) (1) 
b = 2Q(jJ\T^\JJ^ (2) 

c = -n(jj\T ( \ 3) \jjy (3) 
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and write the first-order hyperfine energy as 



W. 
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where 1+ = 1(1 + 1), 



K = F. 



I + , and K + 



J+ = J(J + 1), F+ = F(F 
K(K + 1). 



of a state is given by 



W, 



(2) 
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For the state |1) = |6p 3 / 2 ) of Cs, the second-order hy- 
perfine energy is dominated by the single state \n) = 
|2) = |6pi/2). Moreover, the largest contribution from 
this state is the one associated with the magnetic dipole 
term k = k' = 1. After angular reduction, the second- 
order energy is 



W 



(2) _ 



■h 
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Ei — E2 
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Contributions to the second-order energy from the nu- 
clear quadrupole moment have been evaluated in low- 
est order and found to change the values obtained from 
Eq. (J5J by less than 1%. The fine-structure interval 
E\ — E2 in the denominator is determined as the differ- 

D 3/2 



ence between f D2 the centroid of the 6s 6*1/2 6p Pi 



transition 0,0] and foi the centroid of the 6s 2 Si 
6p 2 P\/2 transition £|, both of which have been measured 



/2 



1/2 

to high precision. One obtains E\ 
1.6609669667(11) x 10 7 MHz. 



E2 — fD2 — Jdi 



III. NUMERICAL ESTIMATES 

Correlation corrections to hyperfine matrix elements in 
alkali-metal atoms are large. Thus, for example, a lowest- 
order Dirac-Hartree-Fock calculation of the hyperfine 
constant 03/2 for the 6^3/2 state of Cs (which is propor- 
tional to the diagonal matrix element (3/2\\T^ ||3/2)) 
leads to a result that is a factor of two smaller than 
the experimental value. One expects (and indeed finds) 
corrections of a similar size for the off-diagonal ma- 
trix element (l/2\\T^ ||3/2> appearing in the numera- 
tor of the expression for the second-order hyperfine en- 
ergy. If we assume that the relative size of the correla- 
tion corrections to the two matrix elements mentioned 
above are the same, then we can determine the ratio 
(l/2||T( 1 )||3/2)/(3/2||r( 1 )||3/2) by means of a lowest- 
order calculation and, using that ratio together with the 
experimental value of a 3 /2, obtain an accurate value for 
(l/2||TW||3/2). That was the strategy used to obtain 



TABLE I: MBPT contributions to the hyperfine constants 
a 1 / 2 and a 3 / 2 (MHz) of the Qp\/2 and 6^3/2 states, re- 
spectively, and to the off-diagonal matrix element gT^ = 
gi (l/2||T (1) ||3/2) (MHz) are presented . The resulting third- 
order hyperfine constants are are compared with experiment. 
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Term 


ai/2 


«3/2 






1st 


160.88 


23.92 


26.97 


), 


2nd 


40.66 


18.84 


-34.15 




Bruck 


84.40 


16.08 


-1.12 


$y 


St Rad 


5.43 


-7.51 


24.85 




Norm 


-1.20 


-0.23 


0.04 




3rd 


88.62 


8.33 


23.77 




Total 
Expt. 


290.17 
291.89 


51.09 
50.29 


16.59 



the values = 401 Hz and W^' = 520 Hz for the 

6p 3 / 2 state of 133 Cs quoted in Ref. pj. (The ratio of 
matrix elements was determined in the non-relativistic 
approximation and did not depend on details of the 6p 
wave function. The non-relativistic approximation is not 
a serious problem, however, since the ratio obtained us- 
ing relativistic Dirac-Hartree-Fock wave functions differs 
from the non-relativistic ratio by less than 5%.) 
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IV. CORRELATION CORRECTIONS 

The estimates made in the previous section de- 
pend on the the assumption that correlation corrections 
to reduced matrix elements of the hyperfine operator 
(jllT^H/) are independent of the total angular momen- 
tum j of the 6pj state. To test that assumption, we car- 
ried out correlated third-order MBPT calculations of the 
three 6pj matrix elements. 

In Table m we give a detailed breakdown of contribu- 
tions to the third-order matrix elements. Formulas for 
the first-, second-, and third-order matrix element are 
given in |5j . We use a modified version of these formulas 
in which: (a) the sum of the second-order matrix element 
and the third-order contribution to the random-phase ap- 
proximation (RPA) is replaced by the exact solution to 
the RPA equations, and (b) all one-electron matrix ele- 
ments in third-order are replaced by their RPA counter- 
parts. These modifications give dipole transition matrix 
elements that are gauge invariant in second- and third- 
order [(J. The third-order hyperfine constants for the 
6pi/2 and 6p 3 /2 states evaluated in this way are within a 
few percent of experiment. Since we use the same method 
to evaluate diagonal and off-diagonal matrix elements, we 
expect the third order off-diagonal matrix element to be 
accurate to a few percent. 

Substituting the third-order off-diagonal matrix ele- 
ment given in Table [Q into Eq. ©, we find W^ 2) = 37.3 



Hz and W \ 



(2) 



48.3 Hz for the 6p 3 / 2 state. Combin- 
ing the second-order corrections with the observed 6p 3 /2 
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hyperfine intervals (MHz) from pj 

W 5 -Wi = 251.0916(20) 
W4-W3 = 201.2871(11) 
W a - W 2 = 151.2247(16), 

we obtain the following values for the hyperfine constants 
(MHz) 

a = 50.28825(23) 
b = -0.4940(17) 
c = 0.00056(7). 

These values agree within error limits with those found 
inQ. 
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